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Abstract 

Strongly nonlinear transport through Diluted Magnetic Semiconductor multiquantum wells oc- 
curs due to the interplay between confinement, Coulomb and exchange interaction. Nonlinear 
effects include the appearance of spin polarized stationary states and self-sustained current oscil- 
lations as possible stable states of the nanostructure, depending on its configuration and control 
parameters such as voltage bias and level splitting due to an external magnetic field. Oscillatory 
regions grow in size with well number and level splitting. A systematic analysis of the charge and 
spin response to voltage and magnetic field switching of II- VI Diluted Magnetic Semiconductor 
multiquantum wells is carried out. The description of stationary and time-periodic spin polar- 
ized states, the transitions between them and the responses to voltage or magnetic field switching 
have great importance due to the potential implementation of spintronic devices based on these 
nanostructures . 

PACS numbers: 85.75.-d, 72.25.Dc, 75.50.Pp, 73.63.Hs 
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I. INTRODUCTION 



Spin injection is one of the aims of spintronicsi thanks to the potential apphcations of 
injectors as spin LED devices, etc. Also quantum state transfer from spin electrons to pho- 
tons by interband transitions is actively investigated^*^'^. One of the most efficient ways of 
spin injection to date^'^ is the use of II- VI dilute magnetic semiconductors (DMSs) that ex- 
hibit the giant Zeeman effect'': they have a conductivity comparable to that of nonmagnetic 
semiconductors, and can boast spin polarizations close to 100% at a small applied magnetic 
field. However, spin-injection experiments in semiconductors enter easily the regime of non- 
linear response^. Different effects could contribute to nonlinear transport and therefore to 
nonlinear spin injection. For example, band bending effects^ in nanostructures give rise to 
a nonlinear current due to the interplay between Coulomb interaction and electron tunnel 
in these confined systems, which have quasi-discrete states. 

Other physical mechanisms inherent to these systems explain their current- voltage char- 
acteristics: for instance a large Zeeman level splitting A in an applied magnetic field B. Re- 
cently, spin transport through DMS diodes^ and multi-quantum well structures (MQWS) 
has been analyzed:^^'^^. These works study nonlinear features of the current (hystere- 
sis, multistability) as a function of the external voltage. Under strong dc voltage bias V, 
electric field domains are formed in MQWS due to the interplay between electron-electron 
interaction and resonant tunneling^. In other sample configurations or, for different doping 
density, there are spin polarized self- sustained current oscillations (SSCOs) and the system 
could behave as a spin oscillator— li^. To tailor the properties of these spin oscillators or 
injectors, it is important to perform a systematic analysis of the transition from stationary 
to time dependent current, in terms of sample configuration, external magnetic field, doping 
density, etc. 

In this paper we analyze the response to voltage {V) or magnetic (B) switching in a 
n-doped dc voltage biased semiconductor MQWS having its first quantum well (QW) doped 
with Mn. Both spin polarized stationary states (SSs) and SSCOs are possible stable states 
of the MQWS for different values of the parameters. Stationary states field profiles consist 
of two electric field domains separated by a domain wall which is a charge monopole^. 
Magnetic field switching requires knowing phase diagrams of the current density J and the 
applied voltage V versus the level splitting A (due to the magnetic field B), and these 
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diagrams are among the results of this paper. The phase diagram of V versus A shows 
regions of stable SSCOs embedded in others of stable SSs. The extension of the SSCO 
regions increases with the number of QWs in the structure. Sudden changes of or i? may 
switch or disconnect SSCOs from an initial stable SS or force the domain wall to change 
its location. The SSCOs are due to periodic triggering of charge dipoles at the Mn-doped 
well and their motion towards the collector^^. Large level splitting induced by B due to the 
exchange interaction provides DMS MQWSs with a new degree of freedom which is absent in 
conventional III-V weakly coupled n-doped semiconductor MQWSs^^. Another important 
difference is that, in the latter, both charge dipoles and monopoles may be triggered at 
the injector (depending on its current-field characteristics: its conductivity if the relation 
between current and field is linear) and both may cause SSCOs^ii^. In these materials 
and for moderate conductivity of the injecting contact, switching the voltage V between 
different SSs involves either upward monopole motion or a dipole-tripole mechanisn>i^ii^ii^. 
For sufficiently large conductivity of the injector, the dipole-tripole mechanism ceases to 
exist and voltage switching involves injection of a charge monopole that moves towards the 
collector until it reaches the QW corresponding to the final stable SS^^. Voltage or magnetic 
switching in II- VI MQWSs always involves dipole nucleation at the Mn-doped QW. 

II. MODEL 

Our sample configuration consists of an n-doped ZnSe/(Zn,Cd,Mn)Se weakly coupled 
MQWS. The spin for the magnetic ion Mn++ is S=5/2 and the exchange interaction between 
the Mn local moments and the conduction band electrons is ferromagnetic in II- VI QWs. 
The energy spectrum corresponding to isolated QWs comprising our weakly coupled 
MQWS has the form Ej + h^k±/{2m*), where m* is the effective mass, k± is the in-plane 
wave vector orthogonal to the growth direction and j = 1, . . . , is the QW subband index. 
In the weak magnetic fields considered here, we disregard Landau-level formation and k±_ is a 
continuous variable^i. Using the virtual crystal and mean field approximations, the exchange 
interaction causes the subband energies to depend on spin in those QWs containing Mn ions: 

= E, T A/2, (1) 
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where 



A{B) = 2J,,NMnSBs (t^b) (2) 

for spin s = ±1/2, and Bs, Jsd, Nun and Teg are the Brillouin function, the exchange 
integral, the density of magnetic impurities and an effective temperature which accounts for 
Mn interactions, respectively^"^. 

We model spin-flip scattering coming from spin-orbit or hyperfine interaction by a phe- 
nomenological scattering time Tgf, which is larger than impurity and phonon scattering 
times: Xgcat < ^sf- Vertical transport in the weakly coupled MQWS is spin-independent se- 
quential tunneling between adjacent QWs, so that when electrons tunnel to an excited state 
they instantaneously relax by phonon scattering to the ground state, with the same spin 
polarization^. Lastly, electron-electron interaction is considered within the Hartree mean 
field approximation. 

The equations governing the model are^»^: the discrete Poisson equation relating the 
two-dimensional spin- up and spin-down electron densities, nf , n~ , respectively, to the aver- 
age electric field —Fi at the ith MQWS period (of lengh /), 

e {Fi - = e {nf + n; - Nn), (3) 

and the rate equations for 



± 



dn^ _ ± _ ± A(n+,n^ ,/i+) 

for i = 1, . . . ,N. For numerical convenience, we have introduced here a smoothed form 
A{nf,n~,fif) of the scattering term used in^^, given by 

n+ 



1 + exp 



7m 



where 7^ is a small smoothing parameter (smaller than 7 = /i/rscatt or than the thermal 
energy) such that, as 7^ — > 0, 



Tnf for fxf > E^^^, 
otherwise, 



±A{nf,n- ^ ±n. + 
which was used by Sanchez et al^ (note that fif > E^^ is equivalent to fif — E^^ > A). 
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In these expressions, /if is the chemical potential at the ith MQWS period and E^^ are 
the spin-dependent subband energies (measured from the bottom of the ith. well): E^i = 
Ej =1= A/2, and E^—Ej for i 7^1. Also, Njj and e are the 2D doping density at the QWs and 
the average permittivity. 

In weakly coupled MQWS, tunneling between adjacent QWs can be treated in leading 
order perturbation theory. Since elastic and inelastic scattering times in the QWs are shorter 
than any other time scale of the problem, we can assume that the electrons in each well are 
in quasiequilibrium between succesive tunneling events and that their temperature is that of 
the lattice. We ignore interwell spin-flip processes, so that currents are carried between wells 
by the two spin subsystems in parallel. Then, as in the case of non-magnetic MQWSs, the 
tunneling current densities across the ith barrier J^i^i can be calculated by the Bardeen 
Transfer Hamiltonian methodr^^^^. See the detailed derivation for non-magnetic MQWSs 
in Ref.-^. The well known resulting expression^^i^^ can be approximated by the formula^: 



xfexpf^^U; 
\ \ m*kBT y 



_eFjJ_ 

l + e '^bt 



(6) 



i = 1, . . . , N — 1, provided that scattering-induced broadening of energy levels is much 
smaller than subband energies and chemical potentials; see Appendix A of Ref.— . The spin- 
dependent "forward tunneling velocity", v^-^^^, is a sum of Lorentzians of width 27, with 
7 = ^/Tscatt (the same value for all subbands, for simplicity), centered at the resonant field 
values i^± = (E±+i-E£)/(e/), 



V 



(/)± 



where % is proportional to the transmission coefficient of the ith barrier?^. For electrons 
with spin ±1/2, the chemical potential iif and the electron densities nf are related by 



± m*kBT , 
nf = n— In 



1 ^ (l^t- Et 
1 + exp 



The voltage bias condition can be written as 



(8) 



E = V- (9) 

i=0 
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Defining Jj-^j+i = Jj^i-^i + Ji^i+i, the total current density J(t) can be calculated as 

1 ^ 

^W = T7-TE'^-m. (10) 

"I" i=0 

Then, time-differencing the Poisson equation, inserting the rate equations for nf in the re- 
sult, and assuming a constant applied voltage {dV/dt = 0), we obtain the following equation 
relating Jj^j+i(t) and J{t) for i = 0, . . . , N: 

dF- 

e^ + J,^,+i = J{t). (11) 

Boundary tunneling currents for i = and i = N are determined by using tunneling 
currents with = n^^^ = Nd/2 (identical normal contacts)^^. 

As initial conditions, we set nf = Nd/2 (normal QWs) and F, = (f)FM, where Fm is a 
reference field corresponding to the first local maximum {Fm, Jm) of the tunneling current 
Ji^i+i{F, nf , n~ , nf^-^^, n^^-^) for nf = nf_i = Nd/2 in a nonmagnetic welU^, and is a 
dimensionless average field defined by 

V 



1{N+1)Fm' 

{N + 1)0 is the dimensionless voltage across the MQWS. 



(12) 



III. RESULTS 

We have considered barrier and QW widths of 10 and 5 nm, respectively, rsf=10^^s 
(normal QW) and lO^^^s (magnetic QW), m*=0.16mo, A^z)=10^°cm-2, e=7.1eo, T=5 K, 
^1=15.76 meV, ^2=61.99 meV, 7=1 meV and 7^=0.1 meVA 

To find the relation between A and B, we used the values g = 2, S = 5/2 and Tcs = T+Tq 
with To = 2 K. The prefactor in can be estimated from Fig. 3 of Ref.— to be 23.26 meV. 
Then, for T = 5 K, we find 

A{B) = 23.26^5/2(5/2.084); (13) 

units of B and A are Tesla and meV, respectively. 

Figure [1] depicts the phase diagram of voltage versus -B-induced level splitting A. We 
observe that the extension of the parameter regions corresponding to SSCOs increases with 
A^, the number of QWs in the structure, and that SSCO regions for a fixed N contain SSCO 
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regions for structures with smaller A^. To the left of the main oscillatory regions there are 
sometimes small oscillatory regions which appear as isolated dots in Fig. [H Whether these 
dots are connected to the main oscillatory regions by extremely thin regions (narrower than 
the discretization error of the code) is something our numerical solution of the model has 
not been able to decide. Even though these connecting regions have not been found and 
therefore they are not shown in the figure, we cannot discard their existence. In fact the 
dots are absorbed by the larger oscillatory regions to their right as increases. 

Figure [2] shows the total current density J(t) as a function of the applied dimensionless 
voltage for = 10. In agreement with Fig. [1], we observe that the width of SSCO 
regions increases with A. For intermediate values of A, the number of oscillatory regions 
first increases and then decreases again when the oscillatory regions merge for larger A. 

The phase diagram shows similar features if we change A at fixed V: there is only one 
finite interval of SSCOs, (A/, A„), cf. Fig. [H In terms of the magnetic field B, the situation 
is somewhat different. The level splitting A is given by a Brillouin function of S, so A can 
only take on values smaller than a saturating value, Aqo- If A^ < Aqo, there may be SSCOs 
for a finite interval of B as shown in Fig. [3] or for an infinite interval if Aqq < A^ (inset of 

Fig. ED. 

The MQWS response to a sudden switching of the voltage or the magnetic field can be 
inferred from Fig. [H Let us increase at fixed A as indicated by the vertical arrow in Fig. [1] 
(from SSCOs to SSs). After a transient, the MQWS settles to the SS, as shown in Figs. ID^a) 
and (b). The periodic generation of a high field domain at the magnetic QW i=l and its 
motion towards the high field region adjacent to the collector yield SSCOs (Panel b). The 
transient corresponds to the nucleation of a last and larger high field domain at the first 
QW and its motion until its trailing domain wall reaches the location corresponding to the 
stable SS (cf. the similar dipole-tripole mechanism ini^iilii^) . Increasing abruptly between 
SSs in different regions of Fig. [T] always involves dipole emission at the magnetic QW, unlike 
the one-well upward domain wall motion possible in conventional III-V MQWS^. 

For fixed 0, a sudden increment of B from a stable SS region to a SSCO region (horizontal 
arrow in Fig. [1]) induces SSCOs, as shown in Fig. Hl^c-d). The transient stage between the 
SS and SSCOs after switching B is due to the formation of a high field domain at the first 
QW which travels towards the collector. After the domain reaches the MQWS end, a new 
high field domain is formed at the first QW and the same situation is periodically repeated. 
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IV. DISCUSSION 



In this paper, we have systematically analyzed the transition from stationary states to self- 
sustained current oscillations through a dilute magnetic semiconductor multi-quantum well 
structure. Switching suddenly a control parameter as the (dimensionless) applied voltage 
or the external magnetic field B may force the system to move between stable oscillatory 
and stationary states through the transition region. Since self-sustained current oscillations 
are caused by triggering high field domains at the magnetic quantum well, we expect our 
results not to change qualitatively with the contact boundary condition. We have used two 
other conditions to check this: 

(i) = n^_^_i = kNd/2 in the tunneling currents for normal contacts, where k is a 
positive constant. 

(ii) The electric field at the injector Fq is calculated by using the Ohm's law: 

d F 

e^ + aF, = J{t), (14) 

instead of the tunneling current formulas^ with known n^. 

The resulting phase diagrams for (i) with different values of k (from 0.5 to 1.5) and for 
(ii) with different values of the contact resistivity 1/cr (from 31.3 to 313 Qm) have the same 
configuration as in Fig. [1] except for small quantitative shifts of the SSCO regions: increasing 
K or l/o" enlarges the oscillatory region. 

This situation is reminiscent of the early theoretical work on the Gunn effect in bulk 
GaAs. Gunn's experiments made it clear that the SSCOs in dc voltage biased n-doped 
GaAs samples are due to the periodic motion of charge dipole waves that appear at the 
cathode and disappear at the anode^^. Theorists soon used a variety of boundary conditions 
at the injecting contact region (cathode) that could produce the required SSCOs mediated by 
charge dipole waves. Among them, Kroemer's contact characteristics (the electron current 
density is a known function of the electric field at the contact: in case (ii), this function is 
linear )21, general models of metal-semiconductor contacts^*^, or contacts with fixed electron 
density but with a notch near the cathode in the doping density profile^^. While the two first 
contact types can give rise to SSCOs mediated by either moving charge dipoles (either high 
or low electric field domains) or monopoles (either charge accumulation or charge depletion 
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layers) depending on parameter values^^*^, a notch in the doping density produces only 
dipoles^Ii^E. 

While it is feasible to list all possible oscillation types in terms of contact parameter values 
(see^^i^), these values cannot be modified once the Gunn diode has been made. Similarly, in 
a conventional III-V weakly coupled n-doped semiconductor superlattice (SL), the boundary 
condition at the injector, the SL configuration and the doping density at the QWs determine 
whether the system exhibits SSCOs mediated by charge dipole or monopole waves^^i^^ or 
multistable static electric field domains; see the review^ and references cited therein. 

In the SL case, there exist partial phase diagrams: (i) doping density vs dc voltage bias 
for fixed boundary condition^, and (ii) injector conductivity vs dc voltage bias for fixed 
doping density (assuming a linear relation between electron current and electric field at 
the contact)^^. However a complete study (which, depending on both doping density and 
injector conductivity, should yield both monopole and dipole SSCOs, as in^-) has not yet 
been carried out. Be this as it may, once the SL has been made and contacted, the stable 
solutions can be selected only by changing the bias and this limits the type of attractors 
present in a particular SL. 

The situation is different in the case of a dilute magnetic semiconductor multi-quantum 
well structure: the magnetic QW plays the role of a "tunable doping density notch" . In prin- 
ciple, any self-oscillations that may appear are due to triggering of dipoles at the magnetic 
QW. However, by changing the external magnetic field we can select either stable stationary 
states or SSCOs as the DMS multi-quantum well response. 

Our results show how to design a device operating a spin injector and a spin oscillator 
by tuning the Zeeman splitting and the parameters determining the sample configuration. 
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FIG. 1: (Color online) Central panel: Phase diagram of dimensionless voltage {N + 1)0 = V/Vo 
{Vq = 0.96 mV) versus A for different A^. Lateral panels: same phase diagram but for specific 
values of = 4, 5, 6, 7, 8, 9, 10. The SS is stable in the white region, whereas SSCOs are stable 
in the shaded (colored) regions. For a given value of A^, the SSCO regions contain those for all 
smaller values of N. 
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FIG. 2: J-V characteristic curves for a MQWS with AT = 10 and level splittings of A = 10, 14 
and 18 meV. 



0.3 




H 5 10 15 20 

Magnetic field B (T) 

FIG. 3: J-B plot for a MQWS with N = 10 and 0=5. Inset: same for = 6, showing indefinite 
persistence of SSCOs for large values of B. 
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FIG. 4: (Color online) J{t) curve and density plot of the electric field F showing the response of 
a = 10 MQWS to the abrupt switches marked in Fig. [TJ (a) and (b): vertical switch from (p=7 
to 9 with A=12 meV. (c) and (d): horizontal switch from B=2 T to 6 T with (j)=5. 
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